Shubnikov-de Haas oscillations of a single layer graphene under dc current bias 
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Shubnikov-de Haas (SdH) oscillations under a dc current bias are experimentally studied on a Hall 
bar sample of single layer graphene. In dc resistance, the bias current shows the common damping 
effect on the SdH oscillations and the effect can be well accounted for by an elevated electron 
temperature that is found to be linearly dependent on the current bias. In differential resistance, 
a novel phase inversion of the SdH oscillations has been observed with increasing dc bias, namely 
we observe the oscillation maxima develop into minima and vice versa. Moreover, it is found that 
the onset biasing current, at which a SdH extremum is about to invert, is linearly dependent on the 
magnetic field of the SdH extrema. These observations are quantitatively explained with the help 
of a general SdH formula. 

PACS numbers: 72.80.Vp, 73.43.Qt, 73.50.Fq 
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The effect of a dc current bias on the nonlinear re- 
sponse of two-dimensional electron systems (2DES) in 
a classically strong magnetic field is a subject of cur- 
rent interest [l|. In conventional 2DES, current bias in- 
duced effects have been widely studied, in the context 
of the breakdown of quantum Hall effect @, and of 
some recently discovered nonlinear phenomena such as 
the Zener-tunneling oscillations [4j and zero differential 
states 0. Nevertheless, similar studies on 2DES with a 
relativistic-like linear energy dispersion, as recently real- 
ized in single layer graphene @, 0], are less reported. 

In this paper, we report on our experimental study on 
the influence of a relatively small dc bias on the magneto- 
transport of a single layer graphene. In the bias regime 
we explored (with current density up to 20 A/m), we 
find the magnetoresistance at lower field (B < 2 T) has 
negligible dependence on dc bias, while the Shubnikov- 
de Haas (SdH) oscillations, occurring at higher fields, are 
obviously damped by increasing bias current. We show 
that the damping of the SdH oscillations can be well ac- 
counted for by an elevated electron temperature that is 
found to be linearly dependent on the bias current. 

Our most important findings, however, are from the 
differential resistance measurements, where a phase in- 
version of the SdH oscillations is observed with increasing 
the bias current. We observe the onset biasing current, 
at which a SdH maxima (minima) is about to invert to 
a minima (maxima), is linearly dependent on the mag- 
netic field of the SdH extrema. These novel observations 
are quantitatively explained by taking into account the 
nonlinear response of the SdH, due to elevated electron 
temperatures by the biasing current. 

Data presented in this paper were measured on a 
lithographically defined Hall bar device of a single layer 
graphene, as shown in Fig. 1(a). The single layer 
graphene was mechanically exfoliated 0, 0] from Kish 
graphite onto degenerately doped silicon substrate with 
a 300- nm thermal oxide SiC>2. The Hall bar pattern was 



defined by electron-beam lithography (EBL) and oxygen 
plasma etching, with PMMA as a resist. Ohmic elec- 
trodes were defined by a second EBL, and by the subse- 
quent 50nm-Pd deposition and lift-off processes. 

Transport measurements were performed on a PPMS 
system (Quantum Design) which can provide a base tem- 
perature of 2 K and a magnetic field up to 14 T. The sam- 
ple was in-situ annealed for an hour to degas the sample 
surface before cooling down. The carrier density of the 
sample was tuned by a gate voltage V g applied to the Sil- 
icon substrate. The differential resistance, r — dV/dl, 
was measured with standard, low frequency (30.9 Hz) 
lock-in technique in the presence of both a small (100 nA) 
ac excitation current and a dc bias current, I^ c ; while the 
dc resistance, R = V/I, was measured by a dc voltage 
meter in the presence of the dc bias current alone. 

The graphene sheet is identified to be single layer by 
the observation of half-integer quantum Hall plateaus, to- 
gether with corresponding minima in magnetoresistance 
R xx , as shown in Fig. lb and Fig. lc. The sample 
mobility is generally higher than 8, 000 cm 2 /Vs within 
the experimental window \V g \ ^ 40 V. The Dirac point, 
V g ~ 1.5 V, is found to be very close to zero gate voltage, 
which indicates the sample is clean. Strong SdH oscilla- 
tions and their very good symmetry about zero magnetic 
field, as shown in Fig. lc, imply high homogeneity of the 
sample. 

In our sample, the most observable effects of a dc bias 
current are on the SdH oscillations. Typical experimental 
traces are shown in Fig. 2, which were measured at T = 
2.0 K and with a fixed gate voltage V g — —40 V. As 
shown in Fig. 2, the magnetoresistance is nearly flat at 
lower field (B < 2 T) and has negligible dependence on 
dc bias, while the SdH oscillations, occurring at higher 
fields, are obviously dampened by dc biasing. 

In dc resistance, the data shown in Fig. 2 (a) resemble 
clearly those of temperature dependence measurements 
shown in Fig. 3 (a), implying an electron heating effect of 




FIG. 1. (color online), (a) A SEM image of the graphene 
Hall bar device. The Hall bar (black color) was defined 
by electron-beam lithography and oxygen plasma etching, 
with six Pd electrodes (white color). The width of the bar 
is about 2.8 /jm, and the distance between voltage con- 
tacts along the same side is about 7 /xm. (b) The lon- 
gitudinal magnetoresistance R xx and Hall conductivity a xy 
against the gate voltage at fixed magnetic field B = 8 T. 
The half- integer quantum Hall plateaus at a xy = ve 2 /h with 
v = 4(N + 1/2), N = ±0, ±1, ±2,. . . , are hallmarks of a sin- 
gle layer graphene. (c) The magnetoresistance measured at a 



fixed gate voltage V g 



-40 V. 



the dc bias commonly observed on the magnetotransport 
of a 2DES. However, except for the amplitude damping, 
the differential resistance shown in Fig. 2 (b) manifests 
a novel feature that the SdH oscillation extrema are in- 
verted with increasing the bias current. Moreover, it is 
found that the onset biasing current (Ii nv ), at which a 
SdH maxima (minima) is about to invert to a minima 
(maxima), is linearly dependent on the magnetic field of 
the SdH extrema (B ex ), with a slope j3 — 4.2 /xA/T, as 
shown in the inset of Fig. 2 (b) . 

. In the regime of SdH oscillations, the magnetoresis- 
tance of a 2DES can be wrote in a general form regardless 
of its energy dispersion[T 
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where A is a constant prefactor, Sf — Ttkp is the area 
enclosed by the Fermi circle, w c = eB/m* is cyclotron 
frequency, r is the lifetime of the carrier, and D(sX) is 




FIG. 2. (color online). The dc bias dependence of mag- 
netoresistance measured at V g = —40 V and T= 2.0 K. (a) 
Traces of dc resistance at selected bias currents, (b) Traces of 
differential resistance at various dc bias currents: from bot- 
tom to top, I dc = 0, 2.0, 4.0, . . . , 40.0 j«A, in 2.0 fj,A steps for 
the lower 21 traces, and Idc = 45.0, 50.0 /iA for the upper 
two traces, respectively. The traces are shifted vertically for 
clarity. In differential resistance, the extrema of the SdH os- 
cillations invert whenever the bias current is sufficiently large; 
the onset traces of the inversion are marked by open and filled 
squares in Fig. 2 (b), for selected SdH maxima and minima, 
respectively. The inset of Fig. 2(b) reveals that the onset bias 
current is linearly dependent on the magnetic field of the SdH 
extrema, with a slope f3 = 4.2 (j,A/T. 



the temperature damping factor 
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In Eq. (1), </>o accounts for the Berry phase of the 
2DES, with <fro — for conventional 2DES and (fro = ir for 
single layer graphene. Due to its linear energy dispersion 
e(fc) = vpfik, the effective mass of a single layer graphene 
is dependent on the carrier density: 



= hkp/vp — (h/vp)y/Tm s 



(3) 



From Eq. (1), the amplitude of the SdH oscillations 
at each extremum (B ex ), neglecting higher harmonics to 
the first order, is given by 



AD(X)exp( ). 
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FIG. 3. (color online), (a) Temperature dependence of the 
magnetoresistance measured at V g = —40 V and zero bias, 
(b) The F(Ae X ,T) vs. T plot for various SdH extrema as 
labeled in the graph. The slope of the solid line corresponds 
to m* = m*, where m* c is calculated by Eq. (3). (c) The 
A SX /D(X) vs. 1/B ex plot that reveals the life time of the 
carrier. The solid line corresponds to r = 53.5 fs. Note the 
vertical intercept of the line at 1/B = indicates A = 2.0. 
The inset reveals that the lifetime is nearly constant in the 
temperature range. 



At sufficient high temperature such that 2ir 2 kBT/huj c 
1, a linear relation on temperature for the quantity 
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follows, which can be used to extract the effective mass 
m*, with the theoretic mass to* calculated by Eq. (3). 

Figure 3 (b) shows the plot F(A ex ,T) vs. T for sev- 
eral SdH extrema of the traces presented in Fig. 3 (a). 
It is evident that data from different extrema collapse 
on the same line at temperature T > 10 K, with a 
slope corresponding to the calculated effective mass to* = 
m * = (fi./t , _F)V 7r71s = 0.0332TO e , where the carrier den- 
sity n s = 3.16 x 10 12 /cto 2 is obtained from the measured 
SdH period and vf = 1.1 X 10 6 m/s adopted from lit- 
eratures 0, [Hj]. This excellent agreement testifies the 
validity of Eq. (1) to describe the SdH oscillations in our 
graphene sample. 

With the effective mass known, Eq. (4) suggests 
that the lifetime can be extracted from the slope of a 
log(A ex /D(X)) vs. l/B ex plot. Such plots for the data 



shown in Fig. 3 (a) are presented in Fig. 3 (c). The in- 
set of Fig. 3 (c) indicates that the lifetime, r ss 54 fs, is 
nearly constant up to T = 50 K, which is consistent with 
the observation of temperature- independent resistance at 
low magnetic fields. 

An interesting result from Fig. 3 (c) is that a prcf- 
actor A = 2.0 is obtained for the graphene sample, 
which is different from that of conventional 2DES where 
A = 4 is theoretically predicted @ and experimentally 
confirmed [Tol Il3j |. There is theoretical implication 11 1 
that A = 2 is intrinsic to the Dirac fermions in graphene, 
however, in addition to our work, more experiments are 
demanding to fully test this point. 

As previously mentioned, the resemblance between the 
data shown in Fig. 2 (a) and those shown in Fig. 3 (a) 
implies an electron heating effect of the dc bias on the 
SdH oscillations. To be more quantitative, the electron 
temperature, T e , can be extracted by fitting the experi- 
mental traces with the SdH formula. In particular, with 
a constant lifetime in the studied regime, Eq. (4) gives 
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where A ex (T e ) and A 



(To) are the amplitude of a SdH 
extremum at B — B ex , measured with or without a 
bias current, respectively, at the same base temperature 
T = To. In Fig. 4 we plot the electron temperature 
T e , extracted via the one parameter fitting to Eq. (6), 
against Idc- The results clearly confirm that T e cx I^c, 
with a slope a = 1.07 K//iA. 

The linear dependence of T e on I^c indicates that the 
energy loss of the electron system, P = Pj OU ie oc I 2 lc oc 
T 2 , implicating an dominant energy dissipation by the 
diffusion of the hot electrons into cold electrodes, rather 
than by the emission of phonons Assuming sim- 

ply the Wiedemann- Franz law, k = CaT ei between the 
thermal and electrical conductivities, the electron tem- 
perature can be estimated from the heat balance be- 
tween the loss by electron diffusion and the joule heat- 
ing V(KVT e ) = Pjouie [141 ]. And further assuming a 
quasi-one-dimensional solution along the Hall bar, the 
electron temperature in the middle of the Hall bar is 
roughly T m « R /{2VZ)I dc , where C = ir 2 k 2 B /{3e 2 ) is 
the Lorenz number, and Rq is the resistance of the Hall 
bar at zero magnetic field. Therefore, we can estimate an 
average electron temperature T e = T m /2 — aide, with 
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Taking the experimental value Rq ~ 700S1, we estimate 
a ~ 1A2K//J.A for the data given in Fig. 2 (a), which 
agrees surprisingly well with the experimental value a ~ 
1.07K/fiA as obtained in Fig. 4. 

Having demonstrated the validity of Eq. (1) for de- 
scribing the SdH oscillations in the single layer graphene, 
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FIG. 4. (color online). The electron temperatures, extracted 
from amplitudes of the SdH oscillations, as function of bias 
current. A linear relation with a slope a = 1.07 K//iA is 
found. The inset shows that, at a given bias current, the 
electron temperature is roughly constant in the regime of SdH 
oscillations. 



and established that the effect of a bias current can be 
taken into account by an effective electron temperature 
T e , now we are ready to focus on the differential resis- 
tance that is given by 
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In 



where R xx is given by Eq. (1) with T = T e {Id c ) 
the experimental regime, we have found that Rq(T) and 
t are near constant with respect to the temperature or 
bias current, it follows 

r xx = Ro [1 + A cos {hSp/eB - tt + <j> )} , (9) 

where higher harmonics of the oscillatory terms have 
been neglected, and the oscillation amplitude is 

dD(X e ) 8T e 



A = A (d(X e ) 



dT e dl dc 



exp( 1 



,(10) 



with X e = 2 r K 2 kBT e /hui c . 

The second term in the bracket of right-hand side of 
Eq. (10) is proportional to Id c , but its sign is negative, 
opposite to the first term, because dD(X e ) / dT e < 0, and 
normally we should have dT e /dId c > 0. As a result, when 
the bias current is sufficiently large, the SdH amplitude 
of the differential resistance can become negative, giving 
rise to a inversion of oscillation extrema. The onset of 
the inversion occurs at 

dD(X e ) dT e 



D(X e )+I dc - 



dT, 81, 



0. 



dc 



(11) 



In our sample, the electron temperature is linear depen- 
dent on bias current, such that the solution of Eq. (11) 
satisfies X e — 1.915, i.e., 

k B T e /frw c = 0.097, (12) 

thus we have the onset current for phase inversion 



which explains well the observed relation as shown in 
Fig. 2 (b). Moreover, substitute the observed coefficient 
a = 1.07 K//zA and the effective mass m* = 0.0332rn e 
into Eq. (13), we get a coefficient (i — 3.67 fiA/T , 
which reasonably agrees with the value (3 = 4.2 /iA/T 
determined from the experimental data. 

From the above analysis, we emphasize that the dc- 
bias-induced inversion of SdH oscillations is unique to 
the differential resistance measurements, unlike that of 
magneto-intersubband oscillations where the inversion 
originates in the dc resistance, as recently discovered in 
double quantum wells [IU [lij]. It is evident that this 
phenomenon in differential resistance is generic in 2DES, 
regardless of their energy dispersion. 

We notice that similar dc-bias-induced inversion of 
SdH oscillations has been observed in conventional 2DES 
of high mobilities fli| . where it is believed that the 
phenomenon cannot be simply described by an eleviated 
electron temperature, rather a nonuniform spectral diffu- 
sion has to be taken into account [l[ . Our work indicates 
that, at least for 2DES in the lower mobility regime, the 
observed phase inversion of SdH oscillations can be well 
accounted for by an eleviated electron temperature. 

In summary, we have studied the influence of a dc bias 
on the magnetoresistance of a single layer graphene. In dc 
resistance, electron temperatures extracted from the am- 
plitude of SdH oscillations manifest an linear dependence 
on the bias current, implicating a dominant heat dissi- 
pation mechanism via electron diffusion. In differential 
resistance, a novel phase inversion of the SdH oscillations 
has been observed, with an onset biasing current that is 
proportional to the magnetic field. 
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T e /a = /IB; with j3 = 0.097He/(k B m*a), (13) 
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